Depinning of an anisotropic interface in random media: The tilt effect 
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We study the tilt dependence of the pinning-depinning transition for an interface described by the 
anisotropic quenched Kardar-Parisi-Zhang equation in 2+1 dimensions, where the two signs of the 
nonlinear terms are different from each other. When the substrate is tilted by m along the positive 
sign direction, the critical force F c (m) depends on m as F c (m) — F c (0) ~ — |m| 1,9( - 1 - ) . The interface 
velocity v near the critical force follows the scaling form v ~ |/| 9 *±(m 2 /|/| e+0 ) with 9 = 0.9(1) 
and (f> — 0.2(1), where / = F — F c (0) and F is the driving force. 

PACS numbers: 68.35.Fx, 05.40. +j, 64.60.Ht 



The pinning-depinning (PD) transition from a pinned 
to a moving state is of interest due to its relevance to 
many physical systems. Typical examples include inter- 
face growth in porous (disordered) media under external 
pressure |l|0|, dynamics of a domain wall under exter- 
nal field jpppfl i dynamics of a charge density wave under 
external field |6|[7j, and vortex motion in superconduc- 
tors under external current ||||]. In the PD transition, 
there exists a critical value, F c , of the external driving 
force F, such that for F < F c the interface (or charge, 
or vortex) is pinned by the disorder, while for F > F Cl 
it moves forward with a constant velocity v, leading to a 
transition across F c . The velocity v plays the role of the 
order parameter and typically behaves as 



v ~ (F - F c 



(1) 



with the velocity exponent 



The interface dynamics in disordered media may be 
described via the Langevin-type continuum equation for 
the interface position h(x,t); 



d t h(x, t) = JC[h] + F + r](x, h). 



(2) 



The first term on the right hand side of Eq. (g) describes 
the configuration dependent force, the second is the ex- 
ternal driving force, and the last, the quenched random 
noise, independent of time, describes the fluctuating force 
due to randomness or impurities in the medium. The ran- 
dom noise is assumed to have the properties, (r)(x, h)) = 
and (r](x,h)r](x',h')) = 2DS d (x - x')S(h - h'), where 
the angular brackets represent the average over different 
realizations and d is the substrate dimension. When 
IC[h] = vV 2 h, the resulting linear equation is called the 
quenched Edwards- Wilkinson equation (QEW) ||l0|| . 

Recently, a couple of stochastic models mimicking the 
interface dynamics in disordered media have been intro- 
duced ]TT| , |l2[ , displaying the PD transition different from 
the QEW universality behavior O] . It was proposed that 



the models are described by the Kardar-Parisi-Zhang 
equation with quenched noise (QKPZ) @, where 



K\h] = v\7 2 h+ -(yhf 



(3) 



The nonlinear term comes from the anisotropic nature 
of disordered medium, thus non-vanishing at the critical 
force F c as shown by Tang et al. jl5| in the context of 
vortex dynamics. When effective pinning force in the ran- 
dom impurity takes the form (A h + s 2 A x ) 2 / 3 /(l + s 2 ) 2/3 , 

1/2 1 /2 

where s = d x h is the local slope and A h and A x 
are the amplitudes of random forces in the h and x di- 
rections, respectively, the nonlinear term is derived by 
expanding the random force in power of the slope s, 
leading to A a (A^ — A x ). Therefore, when > A x 
(Ah < A x ) i.e., when the interface is driven along the 
hard (easy) direction, A is positive (negative), and when 
the medium is isotropic, A^ = A^, the nonlinear term 
vanishes. It has been shown [[j^,[L7| that the interface 
dynamics of the QKPZ equation depends on the sign 
of A, in contrast to the thermal case where the sign is 
irrelevant E3]. 

The critical behavior of the PD transition for the 
QKPZ equation has been thoroughly studied in 1+1 di- 
mensions. For A > 0, the PD transition is continu- 
ous, and the interface at F c is characterized in terms 
of the directed percolation (DP) cluster EsJ spanning 
in the perpendicular-to-the-growth direction [ pdp^ ]. In 
this case, the effective nonlinear coefficient diverges as 
~ (F - Fc)" 1, as F -> F+, and the critical force F c de- 
pends on the substrate-tilt m as 

ll/Kl-a) 



F c (m)-F c (0) ~ -|m| J 



(4) 



where v and a are the correlation length and the rough- 
ness exponent, respectively. These exponents are related 
to one another as |jq ] 

(j> = 2v(l - a) - 9. (5) 

For A < 0, the surface at F c forms a facet with a char- 
acteristic slope s c . The effective nonlinear coefficient is 



insensitive to F. When the substrate-tilt m is smaller 
than s c , the PD transition is discontinuous, and F c is 
independent of m, while the transition is continuous and 
F c increases with m for m > s c . The discontinuous tran- 
sition is caused by the presence of a critical pinning force 
due to both the negative nonlinear term and random 
noise, which is localized. Once this pinning force is over- 
come by increasing external force F, the surface moves 
forward abruptly, yielding the velocity jump at F c . The 
amount of the velocity jump decreases with increasing 
m, and vanishes at the characteristic tilt m c — s c . For 
m > m c , the velocity increases continuously from zero, 
and the PD transition is continuous. Accordingly, the 
characteristic substrate-tilt m c is a multicritical point 

60. 

Since the sign of the nonlinear term is relevant in the 
quenched case, it would be interesting to consider the 
case of the anisotropic QKPZ equation in 2+1 dimen- 
sions, where the signs of the nonlinear terms are alterna- 
tive. Thus, we consider 

d t h = u x a£h + v y d 2 y h + ^ (d x h) 2 + ^{d y hf 

+F + V {x,y,h), (6) 

where A^ > and X y < 0. The anisotropic case is in 
particular interesting due to its application to the vortex 
motion in disordered system |J and the adatom motion 
on step edge in epitaxial surface ]l9[ |. It has been shown 
ijlTjf that for the thermal case, the two nonlinear terms 
cancel each other effectively, thereby the anisotropic 
KPZ equation is reduced to the linear equation, the EW 
equation. For the quenched case, however, the interface 
dynamics in each direction arc different from each other 
and the surface morphology is anisotropic: the surface is 
gently sloping in the positive sign direction (^-direction), 
and is of the shape of a mountain range with steep slope 
in the negative sign direction (y-direction). In spite of 
the facet shape in the negative sign direction, the PD 
transition is continuous due to the critical behavior in 
the positive sign direction. Consequently, one may ex- 
pect that the critical force and the effective nonlinear 
coefficient along the positive sign direction can be de- 
scribed by the scaling theory introduced for A > in 
1+1 dimensions. In this Brief Report, we show, from 
extensive numerical simulations, that this is indeed the 
case and determine the scaling exponents <fi and i>(\ — a) 
independently. 

Direct numerical integration has been carried out using 
standard discretization techniques |20| , ^l[| , in which we 
choose the parameters, v x — v y = 1, A^ = — X y = 1, and 
a temporal increment At = 0.01. The noise is discretized 
as rj{x, y, [h]) 7 where [• • •] means the integer part, and r\ 
is uniformly distributed in [—a/2, a/2] with a = (10) 2 / 3 . 
In order to consider the tilt-dependence, we tilted the 
substrate as h(x, y, 0) = mx along the positive sign 



direction, and used the helicoidal boundary condition, 
h(L + x, y, t) = h(x, y, t) + Lm. We measured the growth 
velocity as a function of the external force F for several 
values of substrate-tilt m, which is shown in Fig. |]. For 
m = 0, we found that 9 = 0.9(1) as in Ref. @. But for 
non-zero m, the exponent 6{m) is generically 1, different 
from its m — value |l5| ]. This picture can be seen in 
Fig. [I], which shows straight lines for large m. 

The critical force F c is estimated as the maximum 
value of F for which all samples of 10 are pinned until 
a large Monte Carlo steps, typically 10 5 Ai. In this way, 
we find F c (0) w 0.51(1). This value is slightly larger 
than that obtained by extrapolating the velocity curve, 
« 0.50. Also, we estimated the critical force F c (m) as a 
function of the substrate-tilt m. The critical force F c (m) 
decreases with increasing substrate-tilt m with the expo- 
nent l/f(l — a) ss 1.9(1) as shown in Fig. g. 

The PD transition is continuous due to the positive 
nonlinear term. To determine the exponent <fi indepen- 
dently, we assume the scaling form for the interface ve- 
locity v(F,m) as 

•H/r* ± (^). (r) 

Here, f = F — F c (0) and the subscript +(— ) denotes 
the positive (negative) / branch. To be consistent with 
Eqs. (0) and (g) and the fact that 6(m) = 1 for m ^ 0, 
the scaling function should behave as, 

*±(a; -» oo) ~ x 6 '^, 
*+(0 + ) = constant, 

and for some positive constant Xq, 

*-N = o, 

^'-(xq) = constant, 

where the prime denotes the derivative with respect to 
x. We find that the velocity data v(F, m) near the tran- 
sition can be collapsed onto a single curve consistent 
with the scaling form. The best collapse is achieved 
with F c (0) w 0.50 and the exponents w 0.9(1) and 
4> w 0.2(1), as shown in Fig. g. Fr om a simple fit to the 
data, we obtain an approximate functional form of the 
scaling function ^±(x) as 

tf + (a;) w A(x + B) 8 K e+ ® (8) 

and, 

*_(x)«/l(x e/(e+<#,) -C), (9) 

with constants A w 0.57(1), B rj 2.1(1), and C = 
x y( e +4>) K 1.5(1). Here we assume that ^ +(x) is ana- 
lytic at x = 0. In Fig. |3[ we also plot Eqs. (g) and (||) 
with dashed lines. The exponents thus obtained satisfy 



2 



the scaling relation, Eq. (|5|), within the errors. 
Alternatively, one may put the scaling form as 

-"•""^'•(sBWff)- < 10 > 

with 

$(x — » oo) ~ a; 6 *, 
= 0, 

$'(— Co) = constant, 
where the positive constant cq is related to xq via 



Co 



c 



With this form, the data can be de- 
scribed by a single scaling function <&(x). The scaling 
plot using the scaling form Eq. ( |To| ) is shown in Fig. ||. 
There, we have shifted the argument x of the scaling 
function <&(x) by Co to make the argument positive. The 
slope of the log-log plot of the scaling function shows a 
crossover from 0.9 to 1.0 with increasing the substrate- 
tilt m, which confirms the fact that 9(m) = 1 for m =/= 0. 

In summary, we have investigated the critical behav- 
iors of the tilted anisotropic QKPZ equation at the PD 
transition. The PD transition is continuous when it is 
tilted along the positive sign direction. It is shown that 
the data can be collapsed onto a single curve, with the 
exponents 9 = 0.9(1) and <fi — 0.2(1). The functional 
form of the scaling function is also numerically deter- 
mined. We have measured, independently, the exponent 
which describes the variation of F c with respect to m, 
and confirmed their consistency. 

This work was supported in part by the Korean Re- 
search Foundation (1999-015-D10070). 
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FIG. 1. The interface velocity as a function of F for various 
m. The curves correspond to the cases of m = 0.0, 0.2, 0.4, 
0.6, 0.8, 0.9, from right to left. 



3 



0.1 



"~6 



0.01 



0.1 



1 



m 



FIG. 2. Log-log plot of F c (0) - F c (m) versus m. The 
dashed line has slope 1.9, drawn for the eye. 
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FIG. 3. Data collapse for the interface velocity using 
Eq. (§) with the exponents 6 = 0.9(1) and = 0.2(1). The 
dashed lines are drawn for the approximate form of the scaling 
functions, Eqs. (0) and M. 
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FIG. 4. Data collapse for the interface velocity using 
Eq. @ with 9 = 0.9(1) and </> = 0.2(1). The dotted (dashed) 
line has slope 0.9 (1.0), drawn for the eye. 
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